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ABSTRACT 

We report measurements of the Sunyaev-ZeFdovich (SZ) effect in three high- 
redshift (0.89 < z < 1.03), X-ray selected galaxy clusters. The observations were 
obtained at 30 GHz during the commissioning period of a new, eight-element 
interferometer - the Sunyaev-ZeFdovich Array (SZA) - built for dedicated SZ 
effect observations. The SZA observations are sensitive to angular scales larger 
than those subtended by the virial radii of the clusters. Assuming isothermality 
and hydrostatic equilibrium for the intracluster medium, and gas-mass fractions 
consistent with those for clusters at moderate redshift, we calculate electron 
temperatures, gas masses, and total cluster masses from the SZ data. The SZ- 
derived masses, integrated approximately to the virial radii, are 1.9 1 0 * 6 x 1O 14 M 0 
for Cl J1415.1+3612, 3.4l°; 6 x 10 14 M o for Cl J1429.0+4241 and 7.21 1 ; 3 x 1O 14 M 0 
for Cl J1226. 9+3332. The SZ-derived quantities are in good agreement with the 
cluster properties derived from X-ray measurements. 
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Subject headings: cosmology: observations — galaxies: individual (NGC5529) 
— clusters: individual (Cl J1415.1+3612, Cl J1429.0+4241, Cl J1226.9+3332) — 
Sunyaev-Zel’dovich Effect — cosmic microwave background — techniques: inter- 
ferometric 


1. Introduction 

Galaxy clusters are the most massive, gravitationally-bound structures in the Universe. 
Over a Hubble time, they form from the rare, high-density peaks in the primordial density 
field on scales of ~ 10 Mpc. As their abundance and evolution are critically dependent on 
cosmology, there is considerable interest in finding clusters at high redshift (z > 1). To date, 
only a few massive clusters at z > 1 have been identified and confirmed, primarily through the 
detection of extended X-ray emission from the hot intracluster medium (ICM) (e.g., Mullis 
et al. 2005; Bremer et al. 2006; Maughan et al. 2006; Stanford et al. 2006). In addition to X- 
ray observations, optical and infrared observations of the cluster member galaxies, and weak 
lensing of background galaxies by the deep cluster potential, are complementary probes of 
high-redshift clusters (e.g., Refregier 2003; Gladders et al. 2006; Stanford et al. 2005; Clowe 
et al. 2006; Wittman et al. 2006). 

Recently, Sunyaev-Zel’dovich effect measurements of galaxy clusters have emerged as a 
powerful probe of cluster physics and cosmology (for a review, see Carlstrom et al. 2002). 
Measurements of the SZ effect have been used to determine cluster properties such as the 
gas and total masses, electron temperatures and scaling relations, as well as to constrain 
Ho and the cosmological distance scale (e.g., Hughes & Birkinshaw 1998; Mason et al. 2001; 
Grego et al. 2001; Reese et al. 2002; McCarthy et al. 2003; Benson et al. 2004; Jones et al. 
2005; Afshordi et al. 2005; LaRoque et al. 2006; Bonamente et al. 2006). Several telescopes 
specifically designed for observations of the SZ effect in high-redshift clusters are currently 
operating or under development, including the Sunyaev-Zel’dovich Array, the Arcminute 
Microkelvin Imager (Kaneko 2006), the Atacama Cosmology Telescope (Fowler 2004) and 
the South Pole Telescope (Ruhl et al. 2004). 

The SZ effect is a spectral distortion of the cosmic microwave background (CMB) radia- 
tion caused by inverse Compton scattering of the CMB photons by electrons in the hot ICM 
(Sunyaev & Zel’dovich 1970, 1972, see also Birkinshaw 1999). The magnitude of the effect 
is proportional to the integrated pressure of the ICM, i.e., the density of electrons along the 
line of sight, weighted by the electron temperature. The SZ flux of a cluster is therefore a 
measure of its total thermal energy. 
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The change in the observed brightness of the CMB caused by the SZ effect is given by 

A J CMB = f(x) f a T n e ^jdl = f(x)y (1) 

Tomb J rn e c 

where Tomb is the cosmic microwave background temperature (2.73 K), <r x is the Thomson 
scattering cross section, k B is Boltzmann’s constant, c is the speed of light, and m e , n e , 
and T e are the electron mass, number density and temperature. Equation (1) defines the 
Compton y-parameter. The frequency dependence of the SZ effect is contained in the term 

m = (z Jii - 4 ) (1 + <5 sz (x, T ')) , (2) 

where x = hu/k B T G MB: h is Planck’s constant, and <5sz is a relativistic correction, for which 
we adopt the Itoh et al. (1998) calculation, valid to fifth order in kBT e /m e c 2 . The SZ effect 
appears as a temperature decrement at frequencies below ~ 218 GHz, and an increment at 
higher frequencies. 

The redshift independence of the SZ effect in both brightness and frequency (the ratio 
A T/T in equation (1) is independent of the distance to the cluster) offers enormous potential 
for finding high-redshift clusters. A cluster catalog resulting from an SZ survey of uniform 
sensitivity will be limited by a cluster mass threshold that is only weakly dependent on 
redshift for z > 0.5 (via the angular diameter distance). To realize the full potential of SZ 
surveys for cosmology will require a thorough understanding of the relationship between SZ 
observables and cluster mass, achievable through detailed cluster observations. 

The SZA is a new, eight-element array of 3.5-meter precision telescopes designed to 
conduct SZ surveys at 30 GHz and detailed cluster observations at 30 GHz and 90 GHz. We 
present observations obtained during the commissioning of the SZA, of three high-redshift, 
X-ray selected clusters: C1J1226. 9+3332 at z = 0.89, Cl J1429. 0+4241 at 2 = 0.92 and 
C1J1415. 1+3612 at z = 1.03. We describe the analysis of the SZ data and compare the 
derived cluster properties with those determined from X-ray observations. We also compare 
the SZA observations of C1J1226. 9+3332 with previous BIMA observations of the SZ effect 
(Joy et al. 2001) in the same cluster. No previous SZ measurements of the other two clusters 
have been reported. 

The paper is organized as follows: In §2 we describe the Sunyaev-Zel’dovich Array, 
followed by a description of the data acquisition, reduction and calibration. The data analysis 
is presented in §3. In §4 we present the results of the analysis and compare with previous 
SZ and X-ray results. Conclusions are given in §5. 
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2. Observations and Data Reduction 
2.1. The Sunyaev-Zel’dovich Array 

The Sunyaev-Zel’dovich Array (SZA) is a new interferometer, equipped with sensitive 
centimeter (26-36 GHz) and millimeter-wavelength (80-115 GHz) receivers, designed specif- 
ically for detecting and imaging the SZ effect in galaxy clusters. In this paper, we report 
results only from the centimeter- wavelength (hereafter 30 GHz) SZA observations. 

An interferometer has sensitivity to angular scales up to the resolution of its shortest 
baseline. To provide a good match to galaxy clusters, which at z > 0.1 subtend several 
arcminutes on the sky, the SZA was designed with small (3.5 meter) antennas, allowing 
close-packed configurations with sensitivity to scales as large as 5' at 30 GHz. The specific 
choice of diameter provides optimum brightness sensitivity to the SZ effect from distant 
clusters, while filtering out the contamination by the intrinsic anisotropy of the CMB at 
larger angular scales. The instantaneous field of view of the SZA is given by the primary 
beam of the 3.5-meter antennas, approximately 12' (FWHM) at 30 GHz. 

Six of the SZA antennas are arranged in a close-packed configuration with separations 
(baselines) ranging from 4.5 — 11.5 m (0.3 — 1.3 kA at 30 GHz), yielding high brightness 
sensitivity, and a resolution of ~ 2', for the detection of the SZ effect in clusters (see Fig. 1). 
Two outer telescopes, one to the north and one to the northeast of the central array, yield 
baselines of up to 65 m (7.7 kA at 30 GHz) to facilitate simultaneous to optimize the surveying 
speed detection of contaminating compact sources at a resolution of ~ 0.3'. This hybrid 
configuration was chosen to optimize the surveying speed of the SZA in an untargeted search. 
Detailed imaging of the SZ effect can be achieved with alternate telescope configurations and 
with SZA observations at 90 GHz. 

The SZA centimeter-wave receivers use cryogenic 26 — 36 GHz high electron mobility 
transistor (HEMT) amplifiers (Pospieszalski et al. 1995), with characteristic receiver temper- 
atures T rx ~ 11 — 20 K. Typical system temperatures, including atmospheric contributions, 
range from ~ 30 — 60 K in the 8-GHz band from 27 to 35 GHz, used for these observations. 
Optical fibers transport the signal (mixed down to intermediate frequencies of 1 — 9 GHz) to 
a wideband hybrid correlator. Sixteen 500 MHz-wide analog bands are further divided into 
seventeen 31.25 MHz channels by 2-bit digital lag correlators, providing spectral discrimina- 
tion as well as high-sensitivity broadband measurements. Additionally, the wide fractional 
bandwidth broadens the Fourier-space (u-v) coverage appreciably, as is evident in Fig. 1. 
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Fig. 1. — Top left: SZA telescope locations. Bottom left: Resulting Fourier-space (u-v) 
coverage for a single track (all baselines). Top middle: u-v coverage for the short-baselines 
only (0 — 2 kA). Bottom middle: u-v coverage for the long baselines only (2 — 8 kA). Top 
right: resulting point spread function (synthesized beam), short-baselines. Bottom right: 
synthesized beam, long baselines. 


2.2. Observations 

Clusters were typically observed for 10 hours about transit, which we refer to as a 
track. Observations of the clusters were interleaved with observations of a strong unresolved 
source every 15 minutes, to monitor variations in the instrumental gain. Cl J 14 15. 1+3612 
was observed for a total of 15 tracks from August 23 to October 8, 2005, using J1331+305 as 
its calibrator. Cl J1429. 0+4241 was observed for a total of 15 tracks between September 30 
and October 28, 2005, also using J1331+305 as its calibrator. Cl J 1226.9+3332 was observed 
for 3 tracks between August 30 and September 4, 2005, using J1310+323 as its calibrator. 

Data were calibrated and excised (flagged) in a data reduction package specific to the 
SZA (see Section 2.3). Data in this paper were taken during the commissioning period of 
the array, during which ~ 44% of the data were flagged. Causes of flagging were: an offline 
antenna (17%); shadowing of one antenna by another (15%); corruption of one 500-MHz 
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band (6%); lack of bracketing calibrator observations (~ 5%), and rare, spurious correlations 
(< 1%). Subsequent improvements to the instrument and observing strategy have reduced 
flagged data to < 5% in good weather conditions. 

In Table 1, we give the pointing center and cluster redshift for each observation (taken 
from Maughan et al. (2006)), along with details of the observations, including the synthesized 
beam sizes (see §3) for both the short and long baselines. We also present the achieved rms 
flux sensitivities for maps made with short and long-baseline data, and the corresponding 
brightness temperature sensitivities for the short-baseline maps. The rms in the short- 
baseline maps is expected to be marginally lower than those of the long baselines (since there 
are 15 short baselines and 13 long baselines). For some of the observations of Cl J1415. 1+3612 
and Cl J1429. 0+4241, however, one of the inner telescopes was offline, resulting in a larger 
number of long baselines than of short baselines. 


2.3. Data Reduction and Calibration 

The SZA data reduction package consists of a suite of MATLAB 1 routines, which consti- 
tute a complete pipeline for flagging, calibrating, and reducing visibility data (see §3) before 
input to higher-level analysis software. In the pipeline, data are converted to physical units 
and corrected for instrumental phase and amplitude variations. Calibration is performed 


1 The Mathworks, Version 7.0.4 (R14), http://www.mathworks.com/products/matlab 


Table 1. Cluster Observations 


Cluster Name z a Pointing Center (J2000) tj nt b T sys c Short Baselines (0-1. 5kA) Long Baselines (2-8kA) 

a 8 (hrs) (K) beam(" X "Z) d +mjy) e B+K) f beam(" X "Z) d <r(mJy) e 


Cl J1415. 1+3612 

1.03 

14^15 m ll s .2 

36° 12 / 04 // 

34.1 

41.9 

115.4x131.2 -61.1 

0.16 

13.5 

15.7x21.4 87.2 

0.16 

Cl J 1429.0+4241 

0.92 

14^29 m 06 s .4 

42°41'10 // 

32.1 

41.7 

109.9x136.9 -60.9 

0.17 

13.6 

15.6x21.3 83.1 

0.16 

Cl J1226. 9+3332 

0.89 

12 /l 26 m 58 s .O 

33°32 / 45 // 

7.6 

42.9 

117.4x125.4 -64.8 

0.38 

32.9 

16.0x21.2 79.2 

0.37 


a Redshifts from Maughan et al. (2006) 

b On source integration time, unflagged data 

c Mean system temperature scaled to above the atmosphere 

d Synthesized beam FWHM and position angle measured from North through East 
e Achieved rms noise in corresponding maps 
Corresponding brightness sensitivity for the short baselines 
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using antenna-based gain solutions; data are flagged if they do not meet designated criteria 
at each reduction step. The pipeline outputs the calibrated, unflagged cluster visibility data. 

The output of the correlator is a complex, dimensionless correlation coefficient — the 
ratio of correlated power to total power for each baseline. To convert these values to physical 
units (K), we first calculate the system temperature of the antennas by comparison with two 
loads of known temperature (a blackbody calibrator load on the telescope, and the CMB). 
System temperature measurements are made at the beginning and end of every scan, where 
scans are defined as the short (typically 5 to 15 minute) observations of either the cluster or 
the calibrator within a track. Data for a scan are written to disk in consecutive 20-second 
integrations, which is the native timescale on which calibration, flagging and Fourier-space 
calculations are performed. Data with high system temperatures (an indication of poor 
weather) , or spurious calibration values (due to readout error of the power sensors) are 
flagged at this stage of the data reduction. 

Absolute calibration is derived from observations of Mars, using fluxes predicted by 
the Rudy model (Rudy 1987). Since Mars is partially resolved on the longest baselines, 
a strong, unresolved source is used to transfer the calibration from the short baselines. 
This absolute flux calibration is performed bimonthly. Measured antenna efficiencies (which 
include aperture efficiencies, atmospheric and system phase noise, correlator efficiencies, and 
other factors) have fluctuated by less than 5% in the period for which the data in this paper 
were taken. The absolute calibration was cross-checked by comparing SZA observations of 
Jupiter to those of WMAP and CBI (Page et al. 2003; Readhead et al. 2004). Based on 
these measurements, we estimate the absolute flux calibration to be better than 10% during 
these observations. 

Bandpasses are measured at the beginning of every track using a flat-spectrum, unre- 
solved radio source with high signal-to-noise in each 31.25-MHz channel. The first and last 
31.25-MHz channel of each 500-MHz band are flagged, as they are corrupted by aliasing in 
the lag correlator. The bandpasses are corrected to remove differential gain (amplitude and 
phase) across the channels within each 500-MHz band, and the absolute flux calibration on 
Mars removes differential gain amplitude across the sixteen 500-MHz bands. This leads to 
a uniform gain calibration across the entire 8 GHz bandwidth. 

The phase from the changing geometric delay as a source is tracked across the sky is 
removed in hardware by inducing a phase shift in the local oscillator signal at each antenna, 
and digitally, in the lag correlator, to account for the bandwidth. Accurate removal of the 
geometric delay, which impacts the dynamic range of the instrument as well as the accuracy 
to which source positions can be measured, requires precise knowledge of each telescope 
location, as well as a measurement of the displacement of the azimuth and elevation axes of 
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the telescopes. These locations have been determined to < 0.1 mm by direct measurements 
and by observations of a set of strong radio, sources across the sky. 

To remove residual instrumental phases, a strong unresolved source is observed every 
15 minutes (over which time the instrumental phase variations are typically much less than 
35°). The interpolated calibrator phase is removed on a per-antenna basis; data are flagged 
if it is not possible to interpolate the phase, either because a calibrator observation is bad or 
because bracketing observations are missing. Data are also flagged if a phase jump of more 
than 35° is observed between calibrator pointings. Observations of the phase calibrator are 
reused to perform a relative amplitude calibration across a given track, to check for any time- 
varying antenna gain. We have seen no indication that the gain on an individual antenna 
changes significantly on timescales less than 12 hours. The gain amplitudes are scaled so 
that the time-average of the gain over a track is the same for all antennas. 


3. Data Analysis 

In the limit where sky curvature is negligible over the instrument’s field of view, the 
response of an interferometer on a single baseline, known as a,, visibility, can be approximated 
by: 

/ +oo p-\~oo 

/ A N (l,m)I{l,m) x exp{— 27rj[ul + vm]}dldm, (3) 

-OO J — OO 

where u and v are the baseline lengths projected onto the sky, l and m are direction cosines 
measured with respect to the (u,v) axes, A^(l, m) is the normalized antenna beam pattern, 
and 1(1, m) is the sky intensity distribution. 

As implied by equation (3), an image of the source intensity multiplied by the antenna 
beam pattern, also known as a dirty map, can be recovered by Fourier transform of the 
visibility data. Note that in addition to modulation by the primary beam, structure in the 
dirty map is convolved with a function that reflects the incomplete Fourier-space sampling 
of a given observation. This filter function is the synthesized beam (see Figure 1), equivalent 
to the point-spread function for the interferometer. The unprocessed maps shown in the first 
four columns of Fig. 2 are examples of dirty maps. A clean map, such as the ones in the 
last column of Fig. 2, is an image from which the full synthesized beam pattern has been 
deconvolved, and the source model reconvolved with a Gaussian fit to the central lobe of the 
synthesized beam. Where a resolution is quoted in this paper, it is this fitted Gaussian that 
is indicated. 

All quantitative analyses described in this paper are the result of simultaneously fitting 
models of the SZ cluster decrement and sources of contaminating emission — both point-like 
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and extended — as detailed below. In all cases, the model is constructed in the sky-plane 
and then multiplied by the primary beam across the field of view. After performing a 
Fourier transform (as given by equation (3)), the resulting model visibilities are compared 
directly to the SZA data. In this way all fitting is done in the Fourier-plane, where the noise 
characteristics and the spatial filtering of the interferometer are well understood; maps are 
used only for examination of the data and to identify cases where contaminating sources are 
present. The two distinct ranges in Fourier-space coverage demonstrated in Figs. 1 and 2 
facilitate removal of sources unresolved by the long-baseline synthesized beam (point sources 
from the perspective of our instrument). Constraints on the point-source flux come primarily 
from baselines including one or both of the outer telescopes, which probe small scales where 
the cluster signal is exponentially damped, while the compact inner array provides sensitivity 
to extended emission. 

The frequency-dependent shape of the primary beam used in the analysis is calculated 
from the Fourier transform of the aperture illumination of each Cassegrain telescope. The 
illumination is modeled as a Gaussian taper, with a central obscuration corresponding to the 
secondary mirror; the validity of this model has been confirmed by holographic measurements 
of the primary mirrors. 


3.1. Point Source Extraction 


When fitting an unresolved radio source, hereafter referred to as a point source, we 
represent it by a delta function, parameterized by the intensity at the band center, / 31 GHz, 
and a spectral index a over our sixteen 500 MHz-wide correlator bands. The point source 


Table 2. Unresolved Radio Sources 


Cluster Field 

# 

RA 

^RA 

DEC 

O'DEC 

d a 

31 GHz Flux ot h 

1.4 GHz flux c 

Ot 



(J2000) 

(s) 

(J2000) 

(") 

0 

(mJy) 

(mJy) 

(1.4/31 GHz) 

Cl J1415. 1+3612 

1 

14 h 15 m ll s .28 

0.12 

+36°12'05".l 

1.2 

0.02 

0.91 ±0.15 

3.15 ±0.21 

0.40 ± 0.06 


2 

14 h 15 m 07 s .43 

0.13 

+36°16'17".0 

1.3 

4.28 

1.01 ±0.17 

1.28 ±0.16 

0.08 ± 0.07 

Cl J 1429.0+4241 

1 

14 fe 28 m 32 s .65 

0.01 

+42°40'20".8 

0.1 

6.26 

26.11 ±0.32 0.42 ±0.17 

43.6 ±2.2 

0.17 ±0.02 


2 

14 h 28"* 46M4 

0.15 

+42°38'11".4 

1.2 

4.77 

1.38 ± 0.21 

4.04 ±0.24 

0.35 ± 0.05 

Cl J1226. 9+3332 

1 

12 h 27 m 18 s .63 

0.08 

+33°32'06".6 

0.7 

4.35 

4.56 ±0.42 1.1 ±1.4 

23.2 ±1.2 

0.53 ±0.03 


a Distance from observation pointing center 
b Spectral index determined from SZA data alone (27-35 GHz) 
integrated FIRST flux at 1.4 GHz 
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intensity at frequency v is then: 

I ps (l,m) = /31GHz x ( W - 0 K m - (4) 

where l r and m' are the coordinates of the point source on the sky. From equations (3) and 
(4) , it can be seen that the visibility amplitude due to a point source is simply its intensity, 
weighted by the normalized primary beam response at the source location. 

The Cl J 1415. 1+3612 and Cl J1429.0+4241 cluster fields were each found to contain two 
point sources, while only one point source was detected in the Cl J1226.9+3332 field (see 
Table 2 and Fig. 2). All detected point sources have counterparts in the 1.4 GHz FIRST 
catalog (White et al. 1997). Spectral indices of sources detected with high signal-to-noise 
can be constrained using SZA data alone (see Table 2, where the limiting signal-to-noise is 
indicated by the poor constraint on the spectral index of the last source listed). The spectral 
indices of weaker sources are not well constrained by our data; for these sources the spectral 
index was fixed to the value determined from the integrated flux at 1.4 GHz and the fitted 
total flux at 31 GHz. The strong source in the Cl J1429. 0+4241 field provides a good test 
of our ability to extract point sources. Following removal of this source from fits to the 
long-baseline data alone, we obtain source flux residuals of less than 3% in both the long 
and short baseline maps. 


3.2. Extended Source Extraction 

Where an extended source of emission is present in the field, its frequency depen- 
dence is modeled as in equation (4), with spatial extent given by an elliptical Gaussian. 
One such object, identified as the edge-on spiral galaxy NGC 5529, was detected in the 
Cl J1415. 1+3612 field. The fitted source parameters are in good agreement with the NVSS 
survey (Condon et al. 1998), which detects a source ~ 2' in extent, at a position angle of 
—64°, with an unresolved minor axis. Comparison of the integrated 31 GHz flux for this 
object (2.55 ± 0.22 mJy) with the integrated flux from the NVSS catalog yields a spectral 
index a iA ghz/ 3 ighz = 0.64 + 0.04, consistent with synchrotron emission. Fig. 2 shows the 
effect of removing both point sources and this extended emission from the Cl J1415. 1+3612 
cluster field. 


3.3. Cluster Parameter Estimation 

In Fig. 3, we present clean maps of the 30 GHz emission of the three cluster fields after 
the removal of the radio sources described above; the SZ decrement is clearly detected for 
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Fig. 2. — Radio source subtraction in our three cluster observations. The rows, from top 
to bottom, depict observations of Cl J1415. 1+3612, Cl J1429.0+4241, and Cl J1226. 9+3332. 
The first column shows the dirty map of the long-baseline data (u-v radii > 2 kA). The 
second column shows the same maps, with fitted point sources removed, as discussed in 
§3.1, resulting in map residuals consistent with noise. The locations of fitted point sources 
are indicated with cross-hairs. The third column shows the dirty map made from the short 
baseline data (u-v radii < 2 kA). The fourth column shows the same maps, following the 
removal of all sources of radio emission. In the top panel (Cl J 1415. 1+3612), we have removed 
the extended emission from NGC 5529, indicated by an ellipse, resolved out in the long- 
baseline map. The resulting cluster decrements reflect the shape of the synthesized beam 
(see §3 and Fig. 1), convolved with the intrinsic cluster profile. The final column shows the 
clean map of the SZ decrement in each cluster. Images are « 21' on a side, nearly twice the 
FWHM of the primary beam. 

each cluster. Also shown in halftone is the corresponding X-ray emission for each field. The 
radio maps are for qualitative comparison only, illustrating the confidence of the detections 
and the alignment of the SZ decrement with the X-ray emission. A quantitative analysis 
of the SZ profile is rendered by fitting a Fourier-transformed model, described next, to the 
visibility data. 
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Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 


Fig. 3. — Sunyaev-Zel’dovich Effect measurements (contours) overlaid on X-ray images 
(halftone) of Cl J1415.1+3612, Cl J1429.0+4241, and Cl J1226.9+3332. The contours are 
set at brightness levels corresponding to integer multiples of ±20 /jK for Cl J1415. 1+3612, 
±20 /rK for Cl J1429. 0+4241 and ±49 /iK for Cl J1226. 9+3332, corresponding to ±1.5 times 
the map rms, respectively (see Table 1); solid contours are negative and dashed contours are 
positive. In each panel, the FWHM of the synthesized beam of the SZ observations is shown 
by the filled ellipse in the bottom left corner. The X-ray images shown here are from the 
XMM EPIC/MOS instrument, with an effective exposure time of 32 ks for Cl J1415. 1+3612, 
66 ks for Cl .J1429. 0+4241, and 45 ks for Cl J1226. 9+3332. The X-ray events are binned in 
3" pixels, and the images have been adaptively smoothed using a threshold of 70 differential 
counts to determine the radii used in the smoothing kernel. 


We model the cluster gas density by a spherical, isothermal /3-model, described by 

/ r 2\ -3/3/2 

n e (r) = n e0 ^ J , (5) 

where the core radius r c and the power law index (3 are shape parameters, and n e0 is the 
central electron number density. The model is a simple parameterization of the gas density 
profile, traditionally used in fitting X-ray (cf. Mohr et al. 1999) and SZ data. Although 
more sophisticated parameterizations have been shown to better reproduce details of the 
density and temperature profiles, data taken with the SZA in its survey configuration lack 
the resolution to constrain models with additional free parameters. Furthermore, this pa- 
rameterization allows a direct comparison with the results of Maughan et al. (2006), who fit 
isothermal /3-models to the X-ray data for the three clusters considered here. 
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The corresponding SZ temperature decrement is given by 


AT {9) = AT 0 


1+ 91 


n 1 _M 
$ 2 \ 2 2 


(6) 


where 9 = r/D A , 9 c = r c / D a, and Z >4 is the angular diameter distance. The temperature 
decrement at zero projected radius, ATq, is related to n e0 by 


AT 0 m e c? 1 1 T(f) 

T C mb f{x)k B a T T e V^D a 9 c T(f - ' 


Best-fit values for the model parameters are determined using a Monte Carlo Markov 
Chain analysis (Bonamente et al. 2004, 2006; LaRoque et al. 2006, and references therein). 
The Markov chains are a sampling of the multi-dimensional likelihood for the model param- 
eters, given the SZ data; the histogram of values in the chain for each parameter is thus 
an estimate of the probability distribution for that parameter, marginalized over the other 
model parameters. In calculations where the fitted parameters ( 9 C , ATo) are combined with 
other variables (such as f gas in the calculation of electron temperatures described in § 4.1), 
each sample in the Markov chain is paired with a random sample of those variables, to obtain 
the probability distribution for the calculated quantity. For any quantity determined from 
the Markov chains, we quote the maximum-likelihood value, with an uncertainty obtained 
by integrating the distribution for that quantity to a fixed probability density, until 68% of 
the probability is enclosed. 

In the analysis described in this and the following sections, two sets of Markov chains 
were generated: one with a weak, uniform prior on 9 C (9 C £ [8".5 — 30"]), and one with a 


Table 3. Resulting Fitted Model Parameters from SZA Observations 


Cluster Name 

8 

0c('T 

scale (kpc/") b 

ARA (") 

ADEC (") 

AT 0 (mK) 

yo (io- 4 ) 

Cl J1415. 1+3612 

2/3 

11.7 

8.06 

— 26.8±i^;| 

-lO.St+x 

-0.42 + 0.10 

0.80 + 0.20 

Cl J 1429. 0+4241 

2/3 

12.4 

7.84 

12.0i? o 2 5 

-4-8+ 4 7 

-0.75 + 0.09 

1.44 + 0.17 

Cl J1 2 26. 9+3332 

2/3 

14.6 

7.77 


16 . 8 +J 

-1.68 + 0.17 

3.28 + 0.33 


a Adopted core radius from X-ray Measurements (Maughan et al. 2006) 
b Assuming H 0 = 70 km/s/Mpc, Qm = 0.3, and Ha = 0.7 
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strong prior (9 C fixed at the value determined from X-ray observations). In both cases, the 
parameter j3 was fixed to 2/3, consistent with X-ray and SZ observations of large cluster 
samples spanning a wide range in redshift (e.g., Mohr et al. 1999; LaRoque et al. 2006) as 
well as X-ray observations of the clusters discussed here (Maughan et al. 2006). 

Table 3 presents cluster model parameters obtained from the Markov chains. The offset 
from the pointing center of the best-fit cluster centroid is given in columns 5 and 6. For 
the isothermal /3-model with fixed /?, although the total SZ fluxes (and therefore cluster gas 
masses; see §4.1) are well constrained by the SZA data alone, the data do not have sufficient 
resolution to separately constrain AT 0 and 9 C . In Table 3 we therefore give constraints on 
the central decrement and central y parameter for fixed (3 = 2/3, and for 9 C fixed to the 
values from Maughan et al. (2006). For the masses presented in §4.1, however, we present 
results derived from fits with only the weak prior on 9 C . 

In Table 4, we compare model parameters from the ~ 8 hours of unflagged SZA obser- 
vations with ~ 42 hours of BIMA data (see Joy et al. 2001) on Cl J 1226.9+3332. Values for 
the cluster locations, central decrements and y-parameters are in excellent agreement. Note 
that the results agree well even though the largest scale probed by BIMA is roughly half 
that probed by the SZA. 


4. Results and Discussion 

In this section, we use the Markov chains of model parameters described in §3.3 to con- 
struct the probability distributions of cluster properties, including the electron temperature, 
gas mass and total mass. We present values integrated to a given overdensity radius, R&( z ), 
defined as the radius at which the mean density of the cluster is related to the critical density 
of the Universe by a fixed density contrast A(z). The density contrast A(z) is assumed to 
scale with redshift like the mean density of a virialized system, as determined from numerical 


Table 4. Comparison of Fit Results for SZA and BIMA 
Observations of Cl J1226. 9+3332 


Instrument 

ARA (") 

ADEC (") 

AT 0 (mK) 

yo (io- 4 ) 

SZA 

BIMA 

1-8±1;| 
0.7 tti 

16.8±?;J 

12.8±t;| 

-1.68 ±0.17 
-1.75 ±0.16 

3.28 ±0.33 
3.41 ±0.31 


Note: Fitted with fixed 

= 2/3 and 0 C 

= 14".6 
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simulations (Bryan &; Norman 1998). 

An estimate of the gas mass in the cluster chn be obtained by multiplying equation (5) 
by p e m p , the mass of the proton weighted by the mean molecular weight of the electrons in 
the gas, and integrating the result to the desired radius: 

[R / r 2\ -3/3/2 

Mg as (R) = fi e m p n e0 J yl + Airr 2 dr. (8) 


The central electron density n e0 is a function of the electron temperature T e and the model 
parameters A Jo, 0 and 6 C , as given by equation (7). 

The total mass of the cluster can be estimated by assuming hydrostatic equilibrium 
(hereafter HSE). For the electron distribution given by equation (5), this approximation 
yields an analytic solution for the total cluster mass contained within a radius R of: 


Mtotal(R) — 


3 k B T e 0 R 3 
Gfimp r 2 c + R 2 ’ 


(9) 


where G is the gravitational constant, pm p is the mean molecular mass of the gas, and r c is 
the core radius, related to 9 C by the angular diameter distance. We adopt a value of 0.3 Z Q 
for the cluster metallicity when calculating both p e and p, consistent with X-ray observations 
of high-redshift clusters (Maughan et al. 2006). 

From equations (7)-(9), we see that if we assume a value for the ratio of the gas mass 
to the total cluster mass, hereafter referred to as the gas-mass fraction, f gas , an estimate of 
electron temperature can be inferred, allowing the masses to be determined without reference 
to an a priori value for T e . We employ this method below to obtain cluster properties from 
the SZ data. For comparison, spectroscopically determined electron temperatures from X- 
ray measurements can be used to estimate the gas masses, total masses, and f gas directly 
from the Markov chains. 


4.1. SZ Derived Cluster Properties 

Here we estimate cluster parameters by assuming a gas-mass fraction, and using the 
SZ data to solve for T e , following Joy et al. (2001) and LaRoque et al. (2003). A previous 
study of a sample of 38 massive clusters obtained a mean of f gas — 0.116 ± 0.005, from 
masses evaluated within a radius of R2500 (distinct from # 2500 ( 2 )) (LaRoque et al. 2006). In 
the calculation of the gas temperature for a single cluster, we therefore adopt a Gaussian 
distribution of f gas with a mean of 0.116 and standard deviation of 0.030, where we have 
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scaled the reported error in the mean by \/37 to approximate the measured distribution of 
.gas-mass fractions. 

We calculate the masses from the Markov chains obtained with a uniform prior on the 
value of 9 C (9 C e [8".5 — 30 "]). For each entry in the chain, we sample the adopted f gas 
distribution and solve for T e from equations (7)-(9), evaluated at # 25 oo for consistency with 
LaRoque et al. ( 2006 ). We use the resulting T e values to obtain estimates for #2500(2) and 
#200(2)- The results, and the gas and total masses calculated at the most likely value for the 
corresponding #a(z) are presented in Table 5. 

For the ACDM cosmology used throughout this paper, the results of numerical simu- 
lations suggest that the density contrast of a cluster at the virial radius — the boundary 
defined by the transition of the dynamical state of the gas from infalling to hydrostatic equi- 
librium — is approximately A ~ 100 at z — 0 (Voigt & Fabian 2006 , and references therein). 
For these clusters, total masses calculated at this density contrast (i.e., at #100(2)) are within 
30 % of those calculated at #200(2)- The regions within the overdensity radii for which we 
quote results thus sample the cluster properties near the cluster core (#2500(2)) and near the 
virial radius (#200(2))- For the three clusters presented in this paper, #200(2) corresponds to 
angular sizes on the order of 1.5' to 2 . 5 / , angular scales well sampled by the short-baseline 
data. 


4.2. Comparison to X-ray values 

The three clusters discussed in this paper have been observed with either the Chandra 
or XMM- Newton observatories, allowing for a direct comparison of SZ derived properties 
with those derived from X-ray data. In the previous section, we calculated masses to radii 
determined self-consistently from the SZ data; for a meaningful comparison with the X-ray 


Table 5. Cluster Masses and ICM Properties Derived from SZ data 



Quantities within i?2soo(z) 


Quantities within i?200(z) 

Cluster Name 0 C T e 

^2500(z) 

Mgas 

M total 

fgas 

^200(z) 

Mg a s Mi 0 i a i 

fgas 

(") (keV) 

(Mpc) 

(10 12 M 0 ) (10 13 m q ) 


(Mpc) 

(10 12 M o ) (10 13 M o ) 


Cl J1415. 1+3612 [8.5,30] 3.T±g| 

0 i7+°-°3 
Ui '-0.04 

q 4+0.7 
d ’ 4 -0.9 

3-i±8:I 

Dll 4+0-024 
U - 114 -0.040 

n 7o+ 0 06 

U *' U -0.09 

29.5lJ;« 18.6l|;| 

nim + 0048 
0 051 

Cl J1429.0+4241 [8.5,30] 5.21J? 

90+0-03 
U.ZA_ 0 Q4 

7 c+1-9 

7.3tl;g 

c-notS 

U.55„o.06 

55.9l« 0 8 6 33.61®;? 

0 1 49+0-052 

Cl J1226. 9+3332 [8.5,30] 8.3lJ;| 

n O1+0.03 
uol -0.03 

18.9±|;2 


o.ioslam? 

1 14+ 0 - 10 
i>14 -0.08 

114.0!]!;° 71.9l? 3 2 3 

0 14>S+ 0 -049 
u.i4o_ 0 040 
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results, however, these should be evaluated at the same physical radii. In all comparisons of 
SZ and X-ray derived masses, we therefore calculate properties to the physical radii for which 
X-ray results have been reported, namely the i?^oo(zj an d R^mlz) va l ues fr° m Maughan et al. 
(2006). 

In Table 6 we present the SZ-derived electron temperature, gas masses, and total masses, 
with corresponding X-ray determinations reproduced from Maughan et al. (2006). We 
present values for two different priors on 6 C (see §3.3): the same weak prior used above, 
and a strong prior on 9 C to facilitate comparison with X-ray derived properties. Note that 
both priors give consistent results, indicating that although SZA data alone cannot break the 
degeneracy of 6 C with other model parameters {J3 and ATo), it can place strong constraints 
on cluster properties resulting from combinations of these parameters. 

The SZ-derived electron temperatures — and therefore total masses — agree within the 
uncertainties with the X-ray values for the two highest-mass clusters. The SZ and X-ray 
derived gas masses for all three clusters are in good agreement. The SZ derived temperature 
for Cl J1415. 1+3612, however, is marginally lower than the X-ray value. We do not believe 


Table 6. Comparison of SZ and X-ray Derived Temperature and Masses 






SZ Derived Quantities 


X-ray Masses 0, 

Cluster Name 

rtX —ray 
R A(z) 

R 

0c b 

T e 

Mgas 

M total 

Te 

Mgas 

M total 


(Mpc) 

n 

(keV) 

(10 12 Mo) 

(10 13 M e ) 

(keV) 

(10 12 m©) 

(10 13 Mq) 

Cl J1415.1+3612 

R2500(z) 

0.23 

11.7 

o 7+0.8 
' 5 *'— 0.7 

6 2+ 1 * 5 
O-^_0.9 

5-6ti:o 

c 7+1-2 
- 0.7 

£ 7+1.2 
D> ' -1.2 





[8.5, 30] 

o 7+O.8 
* —0.8 

5-6iS;! 

5 .Ui ;°3 





R 200(z) 

0.88 

11.7 

3.71“;* 

38.2±1;1 

24.31*3 

s-rtj;? 

38.5±t;| 

38.3tJ 2 4 ° 




[8.5, 30] 

o 7+O.8 
^•'-0.8 

38.9±& s 

24 .otl;i 




Cl J1429.0+4241 

F^2500(z) 

0.26 

12.4 

5 5+ 0 - 7 

io.8il3 


6.21};* 

7 o + l-5 
' — 1.6 

10-511® 




[8.5, 30] 

5.2l“;“ 


8-6±l:l 





I^200(z) 

0.97 

12.4 

5 -Stiio 

62.2l« 1 7 2 

39.711? 

6-2 ±l;g 

42.9t 7 ;* 

44 . 9 li 7 2 * 




[8.5, 30] 

5.2l“;“- 

62.7tfi 7 3 

37.9l®;“ 




Cl J1226. 9+3332 

R 2500(z) 

0.35 

14.6 

S-SRs 

23.7t*;“ 

20.911® 

lo.e!};} 

2i.5l|;| 

25.01U 




[8.5, 30] 

o q+1-8 

23.711? 

20.61®;? 





I^200(z) 

1.29 

14.6 

s-Q±i;| 

130.2±|i;| 

84.7lJ};| 

10 - 6 ±1:1 

119.0l|;| 

102.0li 7 ;| 




[8.5, 30] 

o q+1-8 

130.4l|i;i 

80.311J;? 





Reproduced from Maughan et al. (2006) 
b Prior on core radius, for SZ-derived quantities 
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the discrepancy is significant, given the confidence level for the detection of this cluster, 
uncertainties in the absolute calibration, and possible systematic errors associated with the 
model assumptions. We note, however, that the inconsistency would not be resolved by 
assuming the SZ derived electron temperature is too low. If X-ray spectroscopic temperatures 
are adopted in place of the SZ derived temperatures, the resulting gas masses calculated from 
the SZ data would be marginally inconsistent with the X-ray gas masses. 


5. Conclusion 

We present measurements and analyses of the Sunyaev-Zel’dovich (SZ) effect in three 
high-redshift (0.89 < z < 1.03), X-ray selected galaxy clusters. The observations were 
obtained at 30 GHz during the commissioning period of the Sunyaev-Zel’dovich Array (SZA), 
an eight-element interferometer dedicated to SZ effect observations. The measurements are 
noteworthy in three respects: 1) they extend the redshift range of reported SZ measurements 
to z — 1.03, 2) they extend the low-mass range of reported SZ measurements down to 
~ 2 x 10 14 M©, and 3) with sensitivity to scales as large as ~ 5', the SZA interferometric 
observations provide sensitivity on angular scales larger than the virial radii of the clusters. 

Assuming isothermality and hydrostatic equilibrium for the intracluster medium, and 
gas-mass fractions consistent with those derived from SZ effect and X-ray measurements at 
moderate redshifts, we calculate the electron temperatures, gas masses, and total masses of 
these clusters. The SZ-derived total masses integrated to f? 200 (z) are l-9^cu x 10 14 M© for 
Cl J1415.1+3612, 3.4l°| x 10 14 M e for Cl J1429.0+4241, and 7.2 x 10 14 M o for C1J1226.9 
+3332. These values do not include an overall calibration uncertainty (< 10%) and do 
not account for possible systematic uncertainties associated with the model assumptions. A 
comparison of SZ derived properties to those derived using X-ray data shows good agreement 
between the two methods. 

The data presented here were taken in the SZA survey configuration. Good agreement 
of these measurements with prior X-ray results demonstrates the capability of the SZA for 
probing high-redshift clusters, and bodes well for upcoming SZ surveys. 
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